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We study the phase diagram of the Bose-Hubbard model on the kagome lattice with a broken 
sublattice symmetry. Such a superlattice structure can naturally be created and tuned by changing 
the potential offset of one sublattice in the optical generation of the frustrated lattice. The super¬ 
structure gives rise to a rich quantum phase diagram, which is analyzed by combining Quantum 
Monte Carlo simulations with the Generalized Effective Potential Landau Theory. Mott phases 
with non-integer filling and a characteristic order along stripes are found, which show a transition 
to a superfluid phase with an anisotropic superfluid density. Surprisingly, the direction of the su¬ 
perfluid anisotropy is changing between different symmetry directions as a function of the particle 
number or the hopping strength. Finally, we discuss characteristic signatures of anisotropic phases 
in time-of-flight absorption measurements. 
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Ultracold atoms in optical lattices are prominently 
used to simulate many-body systems in condensed mat¬ 
ter physics [Il-Q . One of the most striking experiments is 
the Mott insulator-superfluid quantum phase transition 
of ultracold bosons in optical lattice built with counter- 
propagating lasers It can be described by the semi¬ 
nal Bose-Hubbard model 0, [1| , where each parameter is 
precisely adjustable in the experiment. With the rapid 
advances in experimental techniques Q the many-body 
physics can now be analyzed on more complex lattice ge¬ 
ometries. On the one hand, the lattice symmetry can be 
reduced by adding additional lasers or tuning their rela¬ 
tive strength, leading to a superlattice structure [iM3, 
which can giv e rise to insulator phases with fractional 
fillings [IJ-llSj. On the other hand, it is also possible to 
enhance the residual entropy of the many-body system 
by using frustrated lattices, which have recently been re¬ 
alized using sophisticated optical techniques |l9l42l| . 


Theoretically, many interesting phases have been pre¬ 
dicted in frustrated lattices such as spin-liquids 2a-l^. 
valence bond solids [i^, string excitation [sSl, ordered 
metals JM 
solids IT 


chiral fractional edge states [32|, and super- 
361. Unfortunately, however, in all these sce¬ 
narios longer range interactions beyond the on-site Bose- 
Hubbard model are assumed, which require dipolar inter- 
actioirs and are experimentally much harder to handle. 
On the other hand, the intriguing interplay between a su¬ 
perlattice and the kagome lattice has never been explored 
before for the on-site Bose-Hubbard model. This is sur¬ 
prising since a superlattice structure can be created natu¬ 
rally in optically generated kagome lattices, and insulat¬ 
ing phases with fractional filling occur without the need 
of longer range interaction as discussed below. We now 
analyze the detailed phase diagram of the Bose-Hubbard 
model on the kagome lattice with a tunable superlattice 
structure. In addition to the fractionally filled ordered 
phases in the quantum phase diagram, the most striking 
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FIG. 1: (a) Potential from Eq. ([T]) of the optical kagome 
superlattice using an enhanced LW laser in the ^-direction 
with 7 = 1.5. Potential along cuts (dashed lines) in the b) 
a 2 -direction and c) ai-direction, respectively, showing the re¬ 
sulting potential offset A/i for sublattice A. 


features are found in the unusual properties of the Bose 
condensed phase, where the anisotropic superfluid den¬ 
sity spontaneously picks a preferred direction depending 
on filling and interaction. This leads to a characteristic 
signature in time-of-flight experiments. 

Let us first consider the optical generation of a kagome 
lattice, which recently has been achieved in experiment 
by using standiirg waves from a long wavelength 1064 nm 
(LW) laser and from a short wavelength 532 nm (SW) 
laser, which are counterpropagating from three 120° di¬ 
rections 2l| . The superposition of the corresponding two 


triangular lattices results in a kagome lattice if the laser 
strengths are exactly equal from all directions. Any slight 
variation of this setup results in a superlattice structure, 
which of course can in turn be used as an additional 
tunable parameter. For example, enhancing the poten¬ 
tial from the LW laser in the ^-direction by a factor 
1 = VeIVq > 1 results in a combined optical potential 


14/Uo = 7^ — 1 -|- 47COs(\/3fc x) cos{ky) + 2 cos{2ky) 
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FIG. 2: Average density fi = {ri\ + ne + n.c)/3 and density 
difference An = uk — (ns + nc)/2 versus /r from QMC with 
T = Ul?iQQ and A = 243 sites at t/U = 0.025 for different 
offsets An/U. 

—2 cos{4ky) — 4 cos{2VSkx) cos{2ky ), (1) 

where k = •\/37r/2ALW in units of the longer wavelength 
Alw = 1064nm. As depicted in Fig. [T] this potential 
leads to an offset of Ay = 4(7 — l)Vb > 0 on one of the 
sublattices A. Note, however, that this offset preserves 
the parity symmetry along the x- and the j/-direction 
and does not increase the unit cell of the kagome lattice, 
which contains three sites. 

Interacting bosons on this lattice can be represented 
by Wannier states, which leads to the well-known Bose- 
Hubbard model for the description of the lowest band in 
second quantized language 

H = -ba^aj)-b - I) 

(bj) * 

^ hi - A/i ^ hj. (2) 

i i^A 

where the nearest-neighbor hopping amplitude t and the 
onsite interaction U are tunable parameters, which de¬ 
pend on the scattering cross-section and the potential 
depth Vb 0. In principle, the potential shift Ay also 
affects the Wannier states and hence other parameters 
in Eq. (ED, but for reasonably small values of Ay these 
higher order corrections can be neglected since they pre¬ 
serve the symmetry of the problem. The chemical poten¬ 
tial y is used to tune the particle number in the grand- 
canonical ensemble. In the following we will use the 
stochastic cluster series expansion algorithm (stI - I^ for 
unbiased quantum Monte Carlo (QMC) simulations of 
this model. In addition the Generalized Effective Po¬ 
tential Landau Theory (GEPLT) provides an analytic 
method to estimate the phase boundaries in an expansion 
of the hopping parameter t/U 

For vanishing hopping amplitude t in the atomic limit, 
the competition between U and Ay can induce several 
incompressible insulating phases. When y is less than 
—Ay, no site is occupied, and the Mott-0 phase is the en¬ 
ergetically favored state. For a larger chemical potential 



FIG. 3: Quantum phase diagram extrapolated to the ther¬ 
modynamic limit obtained from QMC (blue) and multi- 
component GEPLT in second order in t/U (red) at (a) 
Ay/U = 0, (b) Ay/U = 0.25, and (c) Ay/U = 0.5. The 
vertical lines indicate the parameter ranges in Figs. [2] and |4] 
while the horizontal lines are used in Fig. [5] 

—Ay < y < 0, only sublattice A will be occupied with 
one boson per site while the other sites remain empty. 
This phase is therefore 1/3 filled with an order in the 
form of occupied horizontal stripes. Such a 1/3 striped 
density phase (SD) can also occur spontaneously in the 
extended Hubbard model when nearest and next-nearest 
interactions are included [s^. However, longer range in¬ 
teractions are notoriously difficult in optical lattices, so 
that the proposed superlattice is a convenient tool to 
study this phase. For positive values of ^ > 0, the sys¬ 
tem enters the familiar uniform Mott-1 insulator with 
filling factor one. Continuing this analysis for larger y, 
we deduce that SD-n phases with fractional filling factor 
n — 2/3 occur for U{n — 1) — Ay < y < U{n — 1), which 
are separated by Mott-n insulators with integer filling n 
ioT U{n — 1) < y <Un — Ay. 

Both the integer filled Mott-n phases and the fractional 
SD-n phases remain stable for small finite hopping t. As 
shown in Fig. [5] for t = 0.025 U there are plateaus of 
the average density n = {ua + tl-q -b nc)/3 as a function 
of chemical potential, which are characteristic of those 
incompressible phases. In the fractionally filled SD-n 
phases the density difference An = ua — (ub -b nc)/2 
between the sublattices shows plateaus with a value that 
is slightly reduced from unity due to virtual quantum ex¬ 
citations. The plateau states are separated by compress¬ 
ible phases, which are characterized by a finite superfluid 
density, i.e. an off-diagonal order with a spontaneously 
broken U(l) gauge symmetry which will be analyzed in 
more detail below. 

The corresponding phase diagram is mapped out in 
Fig.[3]using large scale QMG simulations. The second or¬ 
der GEPLT approximation is much less demanding and 
agrees quite well with the QMC data, except near the tips 
of the Mott lobes. With increasing offset Ay the frac¬ 
tionally filled SD phases extend over a larger range not 
only in the chemical potential y but also in the hopping 
t. In fact, the SD-I phase for A/r = 0.5 17 is remarkably 
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FIG. 4: Total superfluid density = {p\ + P 2)/2 and su¬ 
perfluid density difference pi = Pi — pi versus Ailing n from 
QMC for 13U = 1200 and N = 243 at t/U = 0.025. Inset: 
schematic illustration of the different mechanisms for positive 
and negative anisotropy parameters. 


stable up to larger values of hopping t than the uniform 
Mott-1 phase. The transitions to the superfluid phase are 
always of second order and can be understood in terms of 
additional condensed particles (holes) on top of the Mott 
states as the chemical potential is increased (decreased). 

One interesting detail in the phase diagram in Fig.[3]is 
the drastic dependence on A/r of the shape of the Mott-0 
phase transition line in the limit of small hopping, which 
changes from linear behavior p(t) = —At for A/i = 0 
to quadratic behavior p(<) = — Ap — 8t^/Afi for large 
Ap. The linear dependence for Ap = 0 can be under¬ 
stood from a competition of chemical potential with the 
kinetic energy, analogously to the quantum melting on 
the triangular lattice 33|. For finite Ap, on the other 
hand, the melting of the Mott-0 phase takes place by ad¬ 
ditional particles on the sublattice A only, which is not 
connected by any flrst-order hopping processes. In the 
limit of small t, the kinetic energy of those particles is 
therefore determined by the second order hopping coef¬ 
ficient t = t^/Ap, which explains the quadratic behavior 
of the phase boundary. The exact shape of the Mott- 
0 transition p = —Ap/2 — t — ^ Ap^ — 4tAp -|- 36t^/2 
can be determined from the single particle energy on the 
super lattice. 


We now turn to the analysis of the order parame¬ 
ter in the superfluid phase. In the QMC simulations 
we determine the superfluid density along the lattice 
vector direction oi = (1,0) using the winding number 
P\ = (^i)/4/3i and analogous for p\ along the lattice 
direction 02 = (l,-\/3)/2 [^. We use a system with 
N = 243 sites and periodic boundary conditions with 
L — 9 unit cells in both the Si and 02 directions, which 
ensures that p® = p| for the perfect kagome lattice. 
Note, that in general the superfluid density is a response 
tensor with four elements p%^, p%y, Py^, Pyy in the x-y- 
coordinate system [d^. Due to reflection symmetry the 
off-diagonal elements pl.y = Py^ = 0 must vanish. The re¬ 
lation to the superfluid densities along the lattice vectors 


is given by pi = pl^ and p| = {pl^ -k 3p^j^)/4. 

In order to analyze a possible anisotropy we consider 
the average superfluid density p\ = {p\ + p|)/2 and the 
difference p^_ = p\ — p| between the two lattice vector 
directions in Fig. |4] as a function of filling n. For finite 
offsets Ap = 0.25 U and Ap = 0.5 U the superfluid den¬ 
sity is indeed anisotropic, but surprisingly also changes 
the preferred direction with increasing filling n. For low 
densities just above the Mott-0 phase the superfluid den¬ 
sity is dominated by virtual hopping processes between 
the A sublattice. As illustrated in the left inset of Fig. 0] 
this virtual hopping process is not possible along the lat¬ 
tice vector tti, which leads to an anisotropic superfluid 
density with p\ < p|. 

When the filling reaches h = 1/3 the superfluid den¬ 
sity drops to zero in the SD-1 phase as expected, but 
then shows the opposite anisotropy pi > p% for n> 1/3, 
which signals a different mechanism: At h = 1/3 the A 
sublattice is completely filled, so that for slightly larger 
densities h > 1/3 excess particles on the B and C sublat¬ 
tices are now responsible for the superfluid density. As 
shown in the right inset of Fig. 01 the B and C sublat¬ 
tices correspond to connected chains along the Si direc¬ 
tion, which are disconnected by occupied A sites. This 
immediately explains why pi > in this case. 

According to this analysis, positive anisotropies pi > 0 
are therefore a hallmark of an off-diagonal U(l) order 
parameter coexisting with a striped density order of a 
filled sublattice A. This situation is reminiscent of a su¬ 
persolid where a stable density order exists on one fllled 
sublattice and excess particles contribute to the super¬ 
fluidity [ 3 ^, with the main difference that in supersolids 
both the U(l) symmetry and the translational symmetry 
are spontaneously broken. Normally supersolid phases 
require longer range interactions beyond on-site, which 
are experimentally difficult to achieve. The creation of 
supersolid-like regions by introducing a superlattice is ex¬ 
perimentally straight-forward, however. Similar to the 
ordinary supersolid, the supersolid-like regions consid¬ 
ered here are also only stable for relatively small hop¬ 
ping, while for larger hopping the ’’ordinary” superfluid 
behavior dominates as we will see below. 

As long as the hopping t is sufficiently small the al¬ 
ternation of anisotropies between Mott and SD phases 
continues as the density is increased due to the same rea¬ 
soning as above. However, this is not the full story since 
for larger hopping t or larger Ailing h the Mott and SD 
phases are not stable, so it is not clear where the different 
regions of positive and negative pl_ are separated. Indeed 
as shown in Fig. 01 the superfluid density does not drop 
to zero for n = 4/3 and Ap/U = 0.25, since the corre¬ 
sponding line is just outside the lobe of the SD-2 phase 
as shown in Fig. [31(b). Also the anisotropy no longer 
changes sign. We And that in the limit of large hopping t 
the overall density becomes irrelevant. The sublattice A 
remains slightly more occupied for all values of n and t. 
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FIG. 5: Superfluid anisotropy parameter I± = {p\ — pi)/(Pi + 
p|) as a function of t/U from QMC for jHU — 2000 and N = 
432 in the parameter range indicated by the horizontal lines 
in Fig. [3l Insets: QMC simulations of the TOF image for 
t/U = 0.0375, PU = 800, N = 243, Ap/[/ = 0.5, and p/U = 
—0.175 (top) and p/U = 0.425 (bottom). 


Since particles on the A sublattice hardly hop in the Oi 
direction, this leads to p\ < p| in the weak coupling limit 
t > U. We call this behavior the ’’ordinary” superfluid, 
in contrast to the supersolid-like regions of the positive 
anisotropy p!_ > 0, which are basically confined between 
the lobes of the SD-n and Mott-n phases. 

To analyze the crossover between different anisotropy 
regions we show the normalized anisotropy parameter 
I± = (Pi — pViKPi + P"iT a function of t/U for dif¬ 
ferent values /i and Ap. in Fig. 0 For small hopping 
the anisotropy parameter I± is positive in supersolid-like 
regions (p = 0) and negative between the Mott-1 and 
SD-2 phase {fi = U — Ap) as discussed above. For larger 
t/U the anisotropy parameters approach small negative 
values in all cases, corresponding to the ordinary super¬ 
fluid. According to the analysis above, the sign-change 
of I± as a function of t coincides with the delocalization 
of the particles on the A sublattice, which start to con¬ 
tribute to the superfluid density in the 02 direction. This 
behavior can be interpreted as a continuous ’’melting” 
of the supersolid-like phase to the ordinary superfluid, 
reminiscent of the melting of the sublattice order in an 
interaction driven supersolid [33 1. 

Anisotropic superfluid densities appear in a variety 
of different systems such as dipolar Bose-Einstein con¬ 
densates with disorder spin-orbit coupled Fermi 

gases [ 4 ^, coupled spin dimer systems ^|, and systems 
with rectangular shape However, an anisotropic su¬ 
perfluidity which is tunable by the isotropic hopping t 
and changes sign when the order on one sublattice melts 
has not been discussed before to our knowledge. 

The observation of the superfluid-Mott transition 
by time-of-flight (TOF) experiments has been pio¬ 
neered many years ago The TOF absorption pic¬ 
ture measures the momentum distribution S{Cl)/N = 


(I ^^1 and turns out to also show a 

clear signature of the anisotropy parameter. To demon¬ 
strate this effect, we used a QMC technique for calculat¬ 
ing the off-diagonal long-range correlation function dur¬ 
ing the loop update [i^ , which allows a direct simulation 
of the TOF absorption signal. As shown in Fig. [S] for 
pi > 0 (upper inset) and for pi < 0 (lower inset) the 
TOF images display a clear signature of the anisotropy, 
which can be used for straight-forward measurements of 
the melting from supersolid-like to ordinary superfluid 
states. 


In conclusion, we analyzed ultracold bosons in a 
kagome superlattice, which can be created and tuned 
by enhancing the long wavelength laser in one direction 
based on recent progress for creating highly frustrated 
lattices 2l|. By using numerical QMC simulations and 
the Generalized Effective Potential Landau Theory, we 
obtained the entire quantum phase diagram including 
Mott phases and fractionally filled charge density phases. 
In the superfluid phase an anisotropic superfluid density 
is found, which changes direction as the overall density or 
the hopping is changed. By tuning the hopping it is possi¬ 
ble to induce a continuous melting from a supersolid-like 
state with a filled sublattice A and positive anisotropy 
parameter p^ > 0 to an ordinary superfluid phase, which 
generically is characterized by a negative anisotropy pa¬ 
rameter p!_ < 0. Both the fractionally filled insulating 
phases and supersolid phases have received much atten¬ 
tion by using models with longer-range interactions 


35l| . Using the superlattice structure proposed in this 


work these phases become experimentally much more ac¬ 
cessible by a simple laser setup instead of using dipolar 
interactions. Moreover, the characteristic signature of 
those effects can be measured in straight-forward TOF 
absorption experiments, without the need of single site 
resolution. In particular, by implementing off-diagonal 
measurements in QMC loop updates, it was possible to 
simulate TOF flight images which show a clear signature 
of the anisotropic superfluid density and the change of 
its direction, when the melting takes place. 
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